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1.  INTRODUCTION 


"Forest  Fire"  is  the  first  and  simplest  of  a  number  of  models  treating  the  initiation  of  detonation  in 
energetic  materials  (Mader  1970,  1979;  Mader  and  Forest  1976;  Lundstrom  1988).  It  is  useful  because 
it  can  be  calibrated  from  readily  available  sensitivity  data.  The  more  recent  models  are  increasingly 
sophisticated  and  capture  more  of  the  micromechanical  phenomenology  involved  in  initiation  of 
heterogeneous  explosives  (Lee  and  Tarver  1980;  Johnson,  Tang  and  Forest  1985;  Tang,  Johnson  and 
Forest  1985).  Often,  however,  Forest  Fire  is  the  only  choice  available  for  a  given  energetic  material.  A 
problem  arises  in  that,  on  the  basis  of  the  assumptions  used  in  its  derivation.  Forest  Fire  is  quite  limited 
in  its  range  of  applicability.  It  can  be  expected  to  adequately  predict  the  distance  of  run  to  detonation  for 
sustained-shock  loading  as  it  is  directly  calibrated  to  do  just  that.  However,  in  the  complex  environments 
associated  with  accidental  initiation,  pulsed-shock  and  finite-rate  compression  (or  ramp-wave)  loading  is 
common  (Staikenberg  et  al.  1989). 

Indeed,  Forest  Fire  is  routineb  employed  in  these  environments,  and  beyond  its  limits(Bowman  et  al. 
1981;  Starkenbeig,  Huang  and  Aibuckle  1984;  Cost  et  al.  1981).  The  consequences  of  this  on  the 
accuracy  of  the  results  achieved  require  clarification,  as  the  literature  does  not  include  assessments  of 
Forest  Fire  performance  for  pulsed  shocks  or  ramp  waves.  The  present  report  is,  in  part,  an  attempt  to 
rectify  this  situation  by  comparing  computational  predictions  with  experimental  results  for  sustained-  and 
pulsed-shock  loading  of  PBX-9404. 

In  the  course  of  pursuing  this  understanding,  it  became  clear  that  several  relatively  simple 
modifications  might  extend  the  applicability  of  Forest  Fire.  These  modifications  include  improved  mixture 
modeling,  use  of  a  more  realistic  reactive  Hugoniot,  and  incorporation  of  surface  area  burning  into  the 
reaction  rate  model.  These  modifications  not  only  improve  the  performance  of  the  model  but  also  serve 
to  illustrate  the  sources  of  the  shortcomings  in  the  original  Forest  Fire  model. 

2.  FOREST  FIRE 

A  brief  discussion  of  the  Forest  Fire  model  is  included  here  to  give  necessary  background.  A  more 
detailed  exposition  has  been  provided  by  Lundstrom  (1988).  Forest  Fire  may  be  considered  to  consist  of 
a  mixture  model,  a  reactive  shock  model,  and  a  reaction  propagation  model. 
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The  two-phase  reacting  mixture  has  been  characterized  by  die  reactant  mass  fraction  (which  varies 
from  1  to  0  as  reaction  proceeds).  However,  the  product  mass  fraction  (which  varies  from  0  to  1)  provides 
a  more  intuitive  description  and  has  been  used  in  the  present  implementation.  The  reactant  and  product 
mass  fractions  are  simply  related.  The  phases  are  assumed  to  be  in  mechanical  and  thermal  equilibrium; 
that  is,  they  are  assumed  to  have  equal  pressures  and  temperatures.  Generally,  the  times  required  to 
achieve  thermal  equilibrium  (via  conduction)  are  significantly  longer  than  those  required  to  achieve 
mechanical  equilibrium  (via  wave  propagation).  Therefore,  the  condition  of  mechanical  equilibrium  can 
be  satisfied  more  often  than  the  condition  of  thermal  equilibrium.  Also  required  are  equations  of  state  for 
each  phase  giving  pressure  and  temperature  as  functions  of  density,  and  specific  internal  energy.  Forest 
Fire  makes  use  of  the  rather  complicated  "HOM"  equations  of  state  for  solid  reactants  and  gaseous 
products. 

The  reactive  shock  is  described  by  a  reactive  Hugoniot,  connecting  the  unreacted  ambient  state  with 
the  Chapman- Jouget  state.  This  Hugoniot  gives  shock  velocity  as  a  linear  function  of  particle  velocity. 
While  the  reactive  Hugoniot  is  a  part  of  the  Forest  Fire  derivation,  and  directly  affects  the  reaction  rate 
law,  it  is  not  generally  used  to  describe  reactive  shock  waves  in  hydrocode  implementations  where  shock 
waves  are  modeled  using  artificial  viscosity.  In  this  case,  the  propagation  rate  law  is  integrated  from  the 
unreacted  state  through  the  distorted  time  scale  associated  with  the  viscous  shock.  The  approach  has  low 
physical  fidelity,  and  results  are  reported  to  be  improved  by  including  the  artificial  viscous  stress  along 
with  the  pressure  in  calculating  the  reaction  rate.  This  has  the  effect  of  increasing  the  amount  or  reaction 
associated  with  the  shock  wave,  which  may  be  lower  than  that  associated  with  the  reactive  Hugoniot. 
Specifically,  this  represents  the  method  that  has  been  implemented  in  the  2DE  code  (Kershner  and  Mader 
1972;  Mader  1979).  During  most  of  the  development  of  a  detonation,  the  amount  of  reaction  associated 
with  the  shock  is  small,  and  the  process  is  dominated  by  downstream  reaction.  The  final  transition  to 
detonation  occurs  rapidly  over  a  short  distance.  Thus,  use  of  this  technique  can  be  expected  to  have  little 
effect  on  predicted  distances  of  run  to  detonation.  In  the  present  implementation,  artificial  viscosity  is  not 
employed,  and  the  reactive  Hugoniot  is  used  explicitly.  It  should  be  emphasized  that,  although  Forest  Fire 
has  no  explicit  ignition  ("hot  spot")  step,  the  reaction  associated  with  the  shock  wave,  whether  determined 
by  integration  or  by  the  reactive  Hugoniot,  plays  this  role. 

The  empirical  reaction  rate  law  is  derived  assuming  that  bulk  reaction  takes  place.  That  is,  the 
reaction  rate  is  proportional  to  the  reactant  mass  fraction  and  "hot  spots,"  and  surface  burning  phenomena 
are  specifically  excluded.  The  pressure  dependence  is  obtained  with  reference  to  the  "Pop  plot,"  which 
represents  commonly  available  sensitivity  data  for  high  explosives  (run  distance  to  detonation  as  a  function 
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of  initial  shock  pressure)  obtained  in  the  "wedge  test"  (Ramsay  and  Popolato  1965).  The  "single  curve 
build-up  hypothesis"  interprets  the  Pop  plot  as  the  path  to  detonation  in  the  pressure-distance  plane  so  that 
the  solutions  obtained  automatically  reproduce  wedge  test  results.  The  derivation  is  simplified  by 
assuming  that  the  pressure  gradient  behind  the  shock  vanishes.  This  condition  is  met  only  when  the 
shock-wave  transitions  to  detonation  as  a  growing  square  wave.  The  derived  Forest  Fire  rate  is  a  function 
of  the  state  variables  immediately  downstream  of  the  shock  wave.  This  is  fitted  as  a  function  of  pressure 
and  applied  at  all  points  downstream  of  the  shock.  It  is  a  peculiarity  of  this  fit  that  the  rate  is  not  defined 
for  pressures  above  the  Chapman-Jouget  value.  The  original  fit  is  a  polynomial  in  pressure  of  as  many 
as  fifteen  terms. 

3.  COMPUTATIONAL  METHODS 

Computations  were  made  using  the  One- Dimensional  Explicit  Shock  (ODES)  code  (Starkenberg  1989). 
The  use  of  a  one-dimensional  code  naturally  limits  the  study  to  consideration  of  one-dimensional 
phenomena.  However,  this  is  quite  sufficient  to  cover  response  to  sustained-shock,  pulsed-shock,  and 
ramp-wave  loading. 

The  most  significant  feature  of  ODES  is  its  explicit  modeling  of  shock  waves.  This  approach  is 
appropriate  whenever  the  flow  between  the  upstream  and  downstream  states  is  quasi-steady.  Generally, 
mass,  momentum,  and  energy  conservation  along  with  overtaking  characteristic  compatibility,  and  other 
appropriate  conditions  are  solved  for  downstream  conditions  with  upstream  conditions  given.  The  shock 
waves  thus  treated  may  be  either  inert  or  reactive.  In  a  reactive  shock,  some  degree  of  reaction  is 
associated  with  the  shock  jump.  This  eliminates  the  requirement  to  numerically  resolve  those  parts  of  the 
reactive  flow  field  which  may  satisfy  the  quasi-steady  condition,  and  be  concentrated  in  a  narrow  region 
downstream  of  the  shock.  In  the  case  of  a  Chapman-Jouget  detonation,  the  entire  reactive  flow  field  is 
unresolved.  However,  it  should  be  noted  that  such  spatial  concentration  does  not  guarantee  that  the 
quasi-steady  condition  is  met.  In  conjunction  with  the  Forest  Fire  model,  the  reactive  Hugoniot  is  used 
to  close  the  system  of  equations  describing  the  shock  wave,  and  is  directly  responsible  for  determining 
the  amount  of  reaction  associated  with  the  shock. 

Shock  waves  may  be  permitted  to  transition  to  reactive  shock  waves,  and  reactive  shock  waves  to 
Chapman-Jouget  or  overdriven  detonations.  In  addition  to  shock  waves,  ODES  also  gives  explicit 
treatment  to  gradient  discontinuities  (such  as  the  heads  and  tails  of  isentropic  waves),  and  allows  them  to 
transition  to  shock  waves.  This  is  useful  for  modeling  finite-rate  compressions. 
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ODES  produces  three  different  plotting  files.  Boundary  plots  which  show  the  positions  of  gradient 
discontinuities,  shock  waves,  reactive  shock  waves,  detonations,  contact  discontinuities  (material 
interfaces),  and  pistons  may  be  obtained.  Profiles  of  density,  internal  energy,  velocity,  mass  fraction,  and 
numerous  other  variables  as  functions  of  position  at  various  times  are  available.  Histories  of  pressure  and 
velocity  at  specified  Lagrangian  stations  may  also  be  produced. 

4.  FOREST  FIRE  PERFORMANCE 

4.1  Sustained-Shock  Stimulus.  Sustained-shock  computations  were  made  using  a  "thick” 
(semi-infinite)  copper  flyer  impacting  a  thick  PBX-9404  sample  to  produce  the  stimulus.  The  velocity 
of  the  flyer  was  varied  to  produce  different  initial  shock  pressures. 

One  such  computation  was  made  with  a  flyer  velocity  of  550  m/s.  The  boundary  plot  in  Figure  1 
shows  the  positions  of  the  shock  wave  propagating  into  the  copper  flyer,  the  interface  between  the  copper 
and  the  PBX-9404,  and  the  reactive  shock  propagating  into  the  PBX-9404  sample  as  functions  of  time. 
Transition  to  detonation  is  evidenced  by  a  change  in  the  slope  of  the  reactive  shock.  Figure  2  shows 
pressure  profiles  at  various  times.  An  approximation  to  a  growing  square  wave  prevails  during  most  of 
the  build-up  process.  This  stands  in  contrast  to  results  predicted  by  other  models  (Lee  and  Tarver  1980; 
Johnson,  Tang,  and  Forest  1985;  Staikenberg  1989)  where  transition  to  detonation  occurs  when  a 
reaction-driven  secondary  compression  overtakes  the  initial  shock  wave.  The  predicted  pressure  histories 
at  several  Lagrangian  stations,  shown  in  Figure  3,  exhibit  poor  agreement  with  experimental  measurements 
(Wackerie  et  al.  1978).  Two  related  aspects  of  the  computed  results  merit  attentioa  First,  the  shock 
pressures  are  too  high,  and  the  shock  propagates  too  rapidly.  Since  shock  waves  are  treated  explicitly 
here,  the  linear  reactive  Hugoniot  determines  the  amount  of  reaction  in  the  shock,  which  may  be  too  great 
Second,  reaction  proceeds  too  promptly  following  shock  passage.  It  appears  that  improvement  might  be 
obtained  by  incorporating  a  model  to  delay  reaction. 

Predicted  distances  of  run  and  times  to  detonation  plotted  as  functions  of  initial  shock  pressure 
are  compared  with  experimental  data  in  Figures  4  and  5,  respectively.  Forest  Fire  was  derived  to  predict 
the  distance  of  run  to  detonation  as  a  function  of  initial  pressure  (or  Pop  plot),  and  it  is  no  surprise  that 
it  does  this  well  (although  the  slope  of  the  computed  Pop  plot  differs  somewhat  from  the  experimental 
value  [Gibbs  and  Popolato  1990]). 
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Boundary  plot  for  the  impact  of  a  thick  copper  flyer  on  a  PBX-9404  sample  at  550  m/s 
showing  the  positions  of  the  shock  wave  propagating  into  the  flyer,  the  interface  between  the 
flyer  and  the  sample,  and  the  reactive  shock  propagating  into  the  sample  as  functions  of  time. 
Transition  to  detonation  is  evidenced  by  a  change  in  slope  of  the  reactive  shock. 
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Figure  2.  Sequence  of  pressure  profiles  showing  detonation  development  in  a  PBX-9404  sample 
following  the  impact  of  a  thick  copper  flyer  at  550  m/s.  The  shock-wave  transitions  to 
detonation  as  a  growing  square  wave. 


6 


Run  to  Detonation  [mm] 


tH 


*fc>4 


* 

*. 


N* 


» 

\ 


Computed 
■  =  Forest  Fire 


x  =  Experiments! 


% 

*  S 


V\ 

*  i 


Initial  Pressure  [GPa] 


ire  4.  Comparison  of  the  experimental  run  to  detonation  versus  initial  shock  pressure  with  Forest  Fire 
predictions  for  PBX-9404.  The  slope  of  the  computed  Pop  plot  differs  somewhat  from  the 
experimental  slope. 
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Figure  5.  Comparison  of  the  experimental  time  to  detonation  versus  initial  shock  pressure  with  Forest 
Fire  predictions  for  PBX-9404.  The  slope  associated  with  the  computed  results  differs 
somewhat  from  that  associated  with  the  experiments. 
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4.2  Pulsed-Shock  Stimulus.  Experimental  data  for  the  critical  impact  of  thin  mylar  flyers  onl9. 1-nun- 
thick  PBX-9404  samples  are  fit  well  by  a  straight  line  in  a  log-log  plot  of  flyer  velocity  versus  flyer 
thickness  (Weingart  et  al.  1984).  This  is  shown  in  Figure  6.  The  critical  velocities  computed  using  Forest 
Fire  correspond  to  transition  to  detonation  occurring  after  19.1  mm  of  run.  They  differ  grossly  from  the 
experimentally  determined  values.  The  predicted  sensitivity  is  much  too  high,  and  tire  variation  appears 
slightly  nonlinear  for  the  thickest  flyers. 

In  the  experimental  study  by  Wackerie,  Rabie,  Ginsberg,  and  Anderson  (1984),  pressure  histories  were 
also  measured  for  a  thin  copper  flyer  impacting  a  PBX-9404  sample  at  the  same  550-m/s  velocity. 
Comparisons  of  the  pressures  measured  at  several  stations  with  results  from  the  corresponding  numerical 
simulation  are  shown  in  Figure  7.  Again,  the  shock  pressures  are  too  high,  and  reaction  proceeds  too 
promptly  following  shock  passage.  This  is  sufficient  to  negate  the  quenching  effects  of  the  rarefactions. 
The  resulting  error  is  biased  toward  safety  when  Forest  Fire  is  used  to  predict  initiation  hazards.  That  is. 
Forest  Fire  predicts  initiations  that  are  not  expected  to  occur. 

4.3  Ramp  Wave  Stimulus.  An  experimental  study  of  the  response  of  an  explosive  to  ramp  waveswas 
conducted  by  Setchell  (1981).  In  this  study,  the  desired  loading  was  produced  by  the  impact  of  a  copper 
flyer  on  a  pyroceram  layer  covering  a  PBX-9404  sample.  The  constitutive  nature  of  pyroceram  causes 
the  impact  shock  wave  to  disperse  into  a  compression  wave  having  a  finite  rise  time.  The  rise  time  can 
be  controlled  by  adjusting  the  thickness  of  the  pyroceram  cover.  In  these  experiments,  evidence  of 
reaction  increased  as  the  rise  time  decreased.  However,  the  effect  on  the  run  or  time  to  detonation  was 
not  determined.  Because  a  constitutive  model  for  pyroceram  was  not  readily  available,  the  experiments 
were  not  directly  simulated.  However,  ramp  waves  were  generated  computationally  by  specifying  the 
pressure  history  at  a  boundary  (piston)  which  is  free  to  move  accordingly. 

Two  such  computations  were  made  to  assess  the  response  of  Forest  Fire  to  ramp  waves.  In  the  first 
of  these,  boundary  conditions  corresponding  to  a  linear  pressure  increase  from  ambient  to  a  final  value 
of  S.O  GPa  over  an  interval  of  1.0  ps  were  applied.  In  the  second,  the  interval  was  cut  in  half.  In  each 
case,  the  boundary  pressure  was  maintained  at  S.O  GPa  after  reaching  that  value.  For  comparison,  another 
computation  was  made  in  which  a  pressure-step  boundary  condition  was  used  to  produce  a  shock  wave 
of  the  same  strength. 
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Figure  6. 


son  of  experimental  critical  fiver  velocity  versus  fiver  thickness  with  Forest  Fire 


redictions  for  initiation  of  PBX-9404  samples  bv  thin  mylar  fivers.  The  predictions  differ 


substantially  from  the  experimental  results. 
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stations  following  the  impact  of  a  0.76-mm-thick  copper  Aver  on  a  PBX-9404  sample  at 
550  m/s.  The  computed  shock  wave  appears  to  incorporate  too  much  reaction  leading  to  rapid 
downstream  reaction. 
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Figure  8  is  a  boundary  plot  comparing  the  wavefront  positions  in  the  three  computations.  The  piston 
trajectories  appear  at  the  left  side  of  the  plot  Distinct  line  styles  are  used  for  isentropic  waves,  shock 
waves  (inert  and  reactive),  and  detonations.  The  initial  slopes  of  the  ramp  wave  trajectories  are  the  same, 
and  correspond  to  the  isentropic  wave  speed  in  the  undisturbed  explosive.  The  isentropic  waves  transition 
to  inert  shock  waves,  reactive  shock  waves,  and  finally  to  detonations.  A  detonation  is  shown  only  when 
reaction  in  the  shock  wave  is  complete.  However,  for  practical  purposes,  transition  to  detonation  is 
identified  as  the  point  at  which  the  slope  of  the  shock  trajectory  assumes  the  value  associated  with  the 
detonation  velocity  (at  which  point  reaction  is  nearly  complete)  rather  than  the  point  at  which  the  line  style 
changes.  The  times  to  detonation  behave  in  the  expected  fashion  qualitatively.  The  shock-wave 
transitions  to  detonation  earliest,  while  the  ramp  waves  transition  first  to  reactive  shock  waves  and  then 
to  detonations  in  order  of  increasing  rise  time.  However,  the  distances  of  run  to  detonation  for  the  shock 
and  fast  ramp  appear  virtually  identical,  while  the  run  distance  for  the  slow  ramp  is  longest  The  effect 
on  the  times  to  detonation  is  stronger  than  the  effect  on  the  run  distances.  Pressure  histories  at  several 
stations  are  shown  in  Figure  9.  These  indicate  that  the  fast  ramp  has  transitioned  to  detonation  at  a 
4.0-mm  run  while  the  shock  wave  has  not 

While  the  computed  results  are,  at  least  in  part,  qualitatively  consistent  with  the  observations  of 
Setchell,  their  predictive  accuracy  cannot  be  assessed.  Indeed,  close  correspondence  to  experimental 
results  should  not  be  anticipated  if  the  effect  of  rise  time  variation  is  on  the  production  of  hot  spots,  a 
mechanism  which  is  not  included  in  Forest  Fire.  Lower  predicted  sensitivity  to  ramp  waves  is  merely  a 
reflection  of  the  fact  that  no  shock  wave  ignites  the  explosive,  and,  as  the  rise  time  increases,  the 
explosive  resides  at  lower  pressures  for  longer  periods.  This  result  should  be  produced  by  any  model  in 
which  reaction  rate  is  a  monotonically  increasing  function  of  pressure. 

5.  MODIFICATIONS 

It  is  of  interest  to  consider  several  modifications  to  the  Forest  Fire  assumptions  and  resulting  model. 
These  include  replacing  thermal  equilibrium  with  adiabatic  reaction,  changing  the  Hugoniot  describing  the 
reactive  shock  wave,  and  replacing  bulk  reaction  with  surface  burning.  Each  of  these  modifications  has 
been  implemented  without  change  in  the  pressure  dependence  of  the  reaction  rate  law. 
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Figure  8.  Comparison  of  predicted  oaths  to  detonation  for  PBX-9404  samples  subjected  to  "slow-ramp," 
"fast-ramp."  and  shock-wave  compressions.  Transition  to  detonation  is  identified  as  the  point 
at  which  the  slope  of  the  shock  trajectory  assumes  the  value  associated  with  the  detonation 
velocity  and  not  where  the  line  style  changes.  The  times  to  detonation  increase  with  increasing 
rise  time  but  the  run  to  detonation  is  shortest  for  the  fast  ramp. 
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Figure  9.  Comparison  of  predicted  pressure  histories  at  four  Lagrangian  stations  for  PBX-9404  subjected 
to  "slow-ramp."  "fast-ramp,"  and  shock-wave  compressions.  The  fast  ramp  produces  the 
highest  pressures. 
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5.1  Mixture  Modeling.  Forest  Fire  has  been  derived  assuming  that  the  reactant  and  product  phases 
are  in  mechanical  and  thermal  equilibrium.  While  the  mechanical  equilibrium  assumption  is  generally 
plausible,  thermal  equilibrium  has  long  been  recognized  as  unrealistic  (Lee  and  Tarver  1980;  Johnson, 
Tang,  and  Forest  1985).  1  ne  conditions  of  mechanical  and  thermal  equilibrium  in  conjunction  with  known 
equations  of  state  for  each  of  the  phases  are  sufficient  to  define  an  equation  of  state  for  the  reacting 
mixture  from  which  mixture  pressure,  and  phase  densities  and  internal  energies  may  be  determined  as 
functions  of  mixture  density  and  internal  energy  and  product  mass  fraction  in  an  iterative  procedure.  This 
may  be  functionally  expressed  as 


P  =P  (p.e>y)  . 


e~eji p,e,y)  , 


Pp=PP(P*-y) . 


and  ep~ep(p,e,y)  , 

where  the  subscripts  r  and  p  refer  to  reactants  and  products  respectively  and  unsubscripted  variables  refer 
to  the  mixture.  Here,  the  product  mass  fraction,  y,  is  used  instead  of  the  original  reactant  mass  fraction, 
1-y,  because  its  progression  from  0  to  1  as  reaction  progresses  provides  a  more  natural  description  of  the 
process. 

An  alternative  has  been  offered  in  conjunction  with  the  JTF  initiation  model  which  substitutes  an 
isentropic  reactant  condition  for  thermal  equilibrium  (Johnson,  Tang,  and  Forest  1985).  This  approach 
allows  the  given  form  of  the  mixture  equation  of  state  to  be  retained  but  suffers  from  the  disadvantages 
of  requiring  time-consuming  iterative  integration  of  the  isentrope  equation  as  well  as  advection  of  the 
initial  state  along  the  isentrope  throughout  the  flow  held.  Further,  this  formulation  is  limited  to  a 
particular  form  of  the  reactant  equation  of  state. 

If  the  partial  differential  equation  describing  conservation  of  energy  of  the  mixture  is  split  between 
the  phases  such  that  thermal  isolation  of  the  reactant  phase  (adiabatic  rear  uon)  is  imposed,  both  the 
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reactant  and  product  internal  energies  may  be  determined  in  the  primary  integration.  For  one-dimensional 
flow,  the  energy  conservation  equations  for  the  mixture,  and  each  phase  are 


Note  that  the  chemical  energy  released  is  treated  as  a  potential  energy  of  the  reactants  so  that  hd  is  a 
positive  constant  and  e  represents  the  thermal  component  of  the  internal  energy.  Since  the  mixture  and 
phase  internal  energies  are  related  by  a  saturation  condition,  only  two  of  the  foregoing  liquations  are 
needed  in  a  computational  implementation. 

In  this  case,  the  mixture  equation  of  state  must  give  the  mixture  pressure  and  phase  densities  as 
functions  of  mixture  density,  phase  internal  energies,  and  product  mass  fraction: 

p  =p  (p.er,ep,y)  , 

p=pr(p,erep.y)  , 

and  pp=pp(p,er,ep,y)  . 

This  can  be  achieved  assuming  only  mechanical  equilibrium  in  conjunction  with  the  phase  equations  of 
state,  which  may  have  any  form. 

The  effect  on  the  pressure  histories  in  the  case  of  a  550-m/s  thick  copper  flyer  is  shown  in  Figure  10. 
Results  obtained  under  the  adiabatic  reaction  assumption  differ  significantly  from  those  obtained  under 
the  thermal  equilibrium  assumption  at  several  stations.  Adiabatic  reaction  produces  a  shorter  run  to 
detonation. 

The  effect  on  the  Pop  plot  is  shown  in  Figure  11.  There  is  a  small  reduction  in  the  predicted  run 
distance.  This  does  not  significantly  alter  the  quality  of  the  agreement  with  the  experiments. 
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Effect  of  the  mixture  model  on  the  predicted  Pop  plot  of  PBX-9404.  The  Pop  plot  is  shifted 
toward  slightly  shorter  runs  with  adiabatic  reaction. 
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The  effect  on  the  response  to  pulsed  shock  waves  is  shown  in  Figure  12.  There  is  a  substantial 
reduction  in  the  predicted  critical  flyer  velocities  which  yields  poorer  agreement  with  experiments. 

Because  of  the  way  Forest  Fire  is  calibrated,  agreement  with  experimental  data  is  better  when  using 
the  original  model.  This  should  not  be  construed  to  imply  that  thermal  equilibrium  is  in  any  way  a 
superior  assumption.  Since  the  mixture  equation  of  state  is  used  in  the  determination  of  the  Forest  Fire 
reaction  rate  coefficients,  strictly  speaking,  these  coefficients  should  be  redetermined  any  time  the  equation 
of  state  is  changed.  However,  the  present  computations  were  made  to  determine  the  effect  of  the  mixture 
model  with  a  fixed  reaction  rate  law. 

5.2  Reactive  Shock-Wave  Modeling.  In  the  present  implementation,  die  reactive  Hugoniot  is  used 
to  close  the  system  of  equations  describing  the  reactive  shock  wave.  It  effectively  determines  the  amount 
of  reaction  in  the  shock.  As  such,  it  constitutes  Forest  Hit’s  ignition  model.  The  reactive  Hugoniot  also 
plays  a  role  in  determining  the  Forest  Fire  reaction  rate.  This  role  has  not  been  considered,  and  the 
reaction  rate  is  treated  as  a  fixed  property  of  the  explosive  as  the  reactive  Hugoniot  is  varied. 

The  use  of  an  inert  shock  in  conjunction  with  the  Forest  Fire  model,  as  suggested  by  Lundstrom 
(1988),  was  not  considered  advantageous.  This  would  require  numerical  resolution  of  the  entire  reactive 
flow  field.  At  least  some  portion  of  the  reactive  region  may  have  a  scale  significantly  smaller  than  that 
of  the  total  flow  field,  and  its  resolution  may  carry  a  significant  computation-time  penalty.  Use  of  a 
reactive  shock  wave  becomes  more  advantageous  as  the  shock  approaches  detonation,  and  includes  greater 
levels  of  reaction. 

The  inert  shock  Hugoniot  for  the  reactants  is  given  by 

U=Ci+sp  . 

Forest  Hre  is  derived  using  a  linear  reactive  Hugoniot  which  connects  the  ambient  state  to  the 
Chapman-Jouget  state. 


where 


U=C;+S,U  , 


Sr—(Dcj—Ci)IUcj  . 
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Figure  12.  Effect  of  the  mixture  model  on  critical  flyer  velocities  for  initiation  of  PBX-9404  samples 
thin  mylar  fivers.  The  critical  velocity  curve  is  shifted  toward  lower  velocities. 
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The  results  computed  using  Forest  Fire  exhibit  a  level  of  reaction  in  low  pressure  shock  waves  which 
is  too  great  to  be  consistent  with  experimental  observations.  More  realistic  representations  are  provided 
by  higher  order  reactive  Hugoniots  which  are  initially  tangent  to  the  inert  shock  Hugoniot,  and  depart 
more  gradually  from  inert  behavior  before  approaching  the  Chapman- Jouget  state.  For  example, 

where 

s»=[ (°cj  ~ci)'ucj  -siI'uCJn~l  • 

The  linear  and  quadratic  (n=2)  reactive  Hugoniots  are  illustrated  in  Figure  13  along  with  the  inert 
Hugoniot 

The  reactive  Hugoniot  appears  explicitly  in  the  Forest  Fire  derivation,  and,  ordinarily,  the  Forest  Fire 
coefficients  would  be  redetermined  any  time  it  is  changed.  Results  presented  here,  comparing  shock 
initiation  computations  for  the  quadratic  Hugoniot  with  those  for  the  linear  Hugoniot  however,  are  for 
the  original  Forest  Fire  rate  law  since  the  intent  is  to  examine  the  effects  of  the  amount  of  reaction 
associated  with  the  shock  only. 

The  results  of  computations  made  for  sustained-shock  loading  are  shown  in  Figures  14  and  15.  The 
pressure  histories  of  Figure  14  indicate  that  use  of  the  quadratic  Hugoniot  reduces  the  shock  pressure  and 
retards  its  progress.  Most  of  this  effect  occurs  during  the  early  portions  of  shock  propagation  and  the  run 
to  detonation  is  essentially  unchanged.  The  Pop  plot,  as  shown  in  Figure  15,  is  only  minimally  altered 
with  generally  shorter  runs  to  detonation. 

In  the  case  of  pulsed-shock  loading,  a  more  substantial  change  in  the  results  is  noted  as  shown  in 
Figure  16.  The  principal  effect  is  to  alter  the  slope  of  the  critical  velocity  curve  such  that  it  is  much 
closer  to  the  experimental  value.  The  general  level  of  agreement  is  still  poor. 

As  previously  noted,  the  sensitivity  of  Forest  Fire  predictions  to  changes  in  the  amount  of  reaction 
in  the  shock  wave  is  relevant  to  the  use  of  artificial  viscosity  in  conjunction  with  the  reaction  rate  law  in 
order  to  describe  reactive  shock  waves.  This  appears  to  produce  less  reaction  than  associated  with  the 
linear  reactive  Hugoniot.  In  the  2DE  code,  the  amount  of  reaction  has  been  increased  by  adding  the 
artificial  viscous  stress  to  the  pressure  for  purposes  of  computing  the  reaction  rate. 


22 


Figure  13.  Unreacted  and  linear  and  Quadratic  reactive  Hueoniots  for  PBX-9404.  Higher  order  Hugoniots 
depart  gradually  from  the  inert  response  at  the  low  velocity  end. 
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Figure  15.  Effect  of  the  reactive  Hueoniot  order  on  the  predicted  Pod  plot  of  PBX-9404.  The  Pop  plot 


is  shifted  toward  slightly  shorter  runs  at  high  initial  pressures  with  the  quadratic  Hugoniot 
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Figure  16.  Effect  of  the  reactive  Hueoniot  order  on  critical  flyer  velocities  for  initiation  of  PBX-9404 
samples  bv  thin  mvlar  fivers.  The  slope  of  the  predicted  critical  velocity  curve  is  close  to  that 
of  the  experimental  curve. 
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5.3  Reaction  Topology  Modeling.  The  Forest  Fire  reaction  rate  model  is  consistent  with  bulk 
reaction.  The  reaction  rate  is  given  by 


y(y*p)m(\-y)F(p)  . 

where  product  mass  fraction,  y,  is  used  to  characterize  reaction  progress  and  F(p)  is  the  Forest  Fire 
pressure  fit.  The  reaction  rate  may  be  generalized  as  the  product  of  function  of  mass  fraction  and  a 
function  of  pressure. 


yO>.p)-S(y)F(p)  • 

Considerable  research  (Howe  et  al.  1976;  Wackerle  and  Anderson  1984)  points  to  the  presence  of 
surface  burning  mechanisms  originating  at  hot  spots  which  delay  the  onset  of  significant  reaction  after 
shock  passage.  Applicability  to  pulsed-shock  loading  problems  may  be  improved  by  adding  surface  area 
dependence  to  the  reaction  rate  model.  The  surface  area  dependence  may  be  identified  with  the 
mass-fraction  dependent  function  while  the  original  pressure  dependence  is  retained.  A  simplification  of 
a  previously  proposed  surface  area  model  (Starkenberg  1989)  may  be  readily  applied  to  the  generalized 
form  of  the  Forest  Fire  rate  law.  In  this  case. 


S(y)  ■  A 


<y/y,)2f3 

[(l-y)/(l-;y,))2/3 

. 


0<y<yt 
y,Zy<  l 


This  model  is  consistent  with  transition  from  outward  "hole  burning”  to  inward  "grain  burning”  at  a 
point  specified  by  a  single  calibration  parameter,  yt,  referred  to  as  the  "transition  mass  fraction."  Small 
values  of  yt  imply  early  transition  to  grain  burning.  This  parameter  can  be  set  with  reference  to  available 
pulsed-shock  response  data.  The  leading  coefficient.  A,  is  determined  by  equating  the  integral  of  the 
surface  area  function  with  respect  to  mass  fraction  to  the  similar  integral  of  the  Forest  Fire  bulk  reaction 
model.  This  gives  A= 5/6,  and  is  intended  to  help  ensure  that  the  sustained-shock  response  of  the  model 
is  not  significantly  altered.  The  surface  area  function  has  been  plotted  for  two  values  of  yr  and  is 
compared  with  the  bulk  reaction  function  in  Figure  17.  The  function  most  closely  conforms  to  bulk 
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Figure  17.  Comparison  of  bulk  and  surface  reaction  functions.  Low  transition  mass  fractions  produce 
closer  conformity  to  the  bulk  reaction  function. 


reaction  for  small  values  of  the  transition  mass  fraction.  Since  the  surface  area  (and,  hence,  the  reaction 
rate)  vanishes  for  unreacted  explosive,  this  model  cannot  be  used  in  conjunction  with  ramp-  v  ave  loading 
unless  an  ignition  model  is  supplied. 

Using  yv  the  surface  burning  model  was  calibrated  to  reproduce  the  experimental  critical  flyer  velocity 
at  a  flyer  thickness  of  0.254  mm.  This  lies  at  the  middle  of  the  range  of  flyer  thicknesses  considered. 
In  conjunction  with  the  linear  reactive  Hugoniot,  the  calibration  gives  y^O^S.  Results  for  the  thinnest 
and  thickest  flyers,  as  well  as  the  0.254-mm  flyer,  are  shown  in  Figure  18.  The  principal  effect  on  the 
critical  velocity  curve  is  to  translate  it  toward  higher  critical  velocities  as  the  transition  mass  fraction 
increases.  Thus,  away  from  the  calibration  point,  agreement  with  experimental  data  is  poor. 

Comparison  of  the  pressure  histories  for  sustained-shock  loading  is  shown  in  Figure  19.  Transition 
to  detonation  occurs  beyond  the  last  station.  The  effect  on  the  Pop  plot,  as  shown  in  Figure  20,  is 
substantial  This  is  primarily  due  to  the  relatively  large  value  of  y(,  which  produces  a  large  delay  in  the 
onset  of  significant  reaction  following  shock  passage.  The  run  to  detonation  associated  with  the  lowest 
initial  pressure  lies  beyond  the  limits  of  the  plot,  and  was  deemed  too  long  to  compute.  Equating  die  area 
under  the  surface  area  function  to  that  under  the  bulk  reaction  function  did  not  successfully  preserve  the 
Pop  plot  for  the  value  yt  of  used.  The  pressure  profiles  in  Figure  21  for  copper  flyer  impact  at  1,000  m/s 
show  that  transition  to  detonation  occurs  by  means  of  an  overtaking  reaction-driven  compression  rather 
than  a  growing  square  wave.  This  is  also  due  to  the  reaction  delay  associated  with  the  surface  area 
mechanism. 

5.4  Combined  Modifications.  The  use  of  different  reactive  Hugoniots  was  seen  to  rotate  the  critical 
velocity  curve  while  the  use  of  the  surface  area  model  in  conjunction  with  different  values  of  transition 
mass  fraction  was  shown  to  translate  it.  It  appears  that  improved  agreement  with  exp.  .imental  data  might 
be  achieved  by  combining  these  two  modifications. 

Another  calibration  of  the  surface  area  model  using  the  quadratic  reactive  Hugoniot  givesy/=0.0794. 
This  produces  substantially  better  results  for  pulsed-shock  response,  as  shown  in  Figure  22.  It  appears 
that  a  still  higher  order  Hugoniot  and  smaller  transition  mass  fraction  would  provide  further  improvement. 
It  is  notable  that  the  best  results  are  achieved  with  a  very  small  value  of  the  transition  mass  fraction 
corresponding  to  a  very  early  transition  to  grain  burning. 
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Figure  18.  Effect  of  reaction  topology  on  critical  flyer  velocities  for  initiation  of  PBX-9404  samples 

thin  mylar  flyers.  The  transition  mass  fraction  has  been  calibrated  at  >>,=0.55  to  produce 
agreement  with  the  experiments  for  a  flyer  thickness  of  0.254  mm.  The  slope  of  the  predicted 
critical  velocity  curve  is  not  close  to  that  of  the  experimental  curve. 
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Figure  20.  Effect  of  reaction  topology  on  the  predicted  Pop  plot  of  PBX-9404.  The  slope  of  the  Pop  plot 
is  substantially  altered  using  surface  burning  with  y,=0.55. 
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Figure  21.  Sequence  of  pressure  profiles  showing  detonation  development  in  a  PBX-9404  sample 
following  the  impact  of  a  thick  copper  flyer  at  1.000  m/s  using  surface  burning  with  yf=0.55. 
The  shock- wave  transitions  to  detonation  on  being  overtaken  by  a  reaction-driven  compression 
wave,  and  not  as  a  growing  square  wave. 


33 


Flyer  Velocity  [mnv^ 


Experimental 
*  =  weingart  et  al. 


Computed 

■  =  Linear  Hugoniot  -  yt=0.55 
7  •  =  Quadratic  Hugoniot  -  yt=0.0794 

-h-t— i  "i  n  ""T  i  i  i  i  |  e  1,1  i"'  i  t  "I’n  i  ti  | 

her9  ter1  1  cf 

Flyer  Thickness  [mm] 


Figure  22.  Effect  of  combining  a  quadratic  reactive  Hugoniot  with  surface  burning  on  critical 

velocities  for  initiation  of  PBX-9404  samples  by  thin  mylar  fivers.  The  transition  mass 
fraction  has  been  calibrated  at  yt  =0.0794  to  produce  agreement  with  the  experiments  for  a 
flyer  thickness  of  0.254  mm.  The  slope  of  the  predicted  critical  velocity  curve  is  closer  to  that 
of  the  experimental  curve. 


The  effect  on  the  Pop  plot  is  shown  in  Figure  23.  The  improved  agreement  with  experiment  is 
attributable  to  closer  conformity  to  the  bulk  reaction  model  when  using  small  values  of  use  of  yr  The  fact 
that  agreement  with  experiment  is  improved  with  respect  to  that  achieved  using  the  original  Forest  Fire 
model  may  be  regarded  as  fortuitous. 

The  effect  on  the  prediction  of  gauge  data  for  sustained-shock  loading  is  shown  in  Figure  24.  The 
comparison  is  still  not  favorable.  The  amount  of  reaction  in  the  shock  wave  has  been  reduced  but  the 
reaction  delay  following  shock  passage  appears  too  great.  The  presence  of  a  strong  compression  wave 
overtaking  the  initial  shock  is  evident. 

6.  SUMMARY 

Forest  Fire  is  limited  in  its  applicability  beyond  prediction  of  run  to  detonation  for  sustained-shock 
loading.  The  reactive  flow  field  is  inadequately  represented,  even  in  this  case.  Forest  Fire’s  prediction 
of  response  to  pulsed-shock  loading  is  grossly  in  error.  Some  aspects  of  the  predicted  response  to 
ramp-wave  loading  are  qualitatively  correct  but  there  is  no  reason  to  believe  that  Forest  Fire  can  be 
accurately  applied  in  this  case. 

The  assumption  of  adiabatic  reaction  (or  isentropic  reactants)  is  more  physically  realistic  than  that  of 
thermal  equilibrium.  The  effect  of  this  assumption  on  flow  field  prediction  may  be  substantial,  especially 
in  the  case  of  pulsed-shock  loading.  Implementation  of  adiabatic  reaction  through  the  energy  conservation 
equations  is  an  economical  computational  approach. 

Reaction  occurring  within  the  shock  wave  constitutes  Forest  Fire’s  ignition  model.  The  amount  of 
such  reaction  is  determined  either  by  integration  of  the  reaction  rate  law  through  a  viscous  shock  or  by 
use  of  a  reactive  Hugoniot  in  conjunction  with  an  explicit  shock  model.  The  influence  of  the  ignition 
reaction  on  the  response  to  pulsed-shock  loading  is  significant. 

The  inclusion  of  a  surface  burning  model  also  has  a  profound  effect  on  the  computational  results. 
When  combined  with  the  quadratic  reactive  Hugoniot,  it  provides  a  good  representation  of  the  responses 
to  a  pulsed-shock  stimulus.  The  best  results  are  achieved  with  calibration  corresponding  to  a  very  early 
transition  to  grain  burning. 
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The  effects  of  modifications  to  Forest  Fire  are  most  significant  with  respect  to  the  predicted  response 
to  pulsed  shock  waves.  Each  of  the  modifications  considered  produced  large  changes  in  the  critical  impact 
velocities  for  thin  flyers.  The  effects  on  the  Pop  plot,  except  in  the  case  of  surface  burning,  are  not 
significant. 
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23.  Effect  of  combining  a  quadratic  reactive  Hugoniot  with  surface  burning  on  the  predicted  Pod 
plot  of  PBX-9404.  The  predicted  Pop  plot  conforms  more  closely  to  the  experimental  plot 
using  the  quadratic  Hugoniot  and  surface  burning  with  >>,=0.0794  than  it  does  using  the 
original  Forest  Fire  assumptions. 
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Figure  24.  Comparison  of  experimental  pressure  histories  with  Forest  Fire  predictions  at  four  I  j* 
stations  following  the  impact  of  a  thick  copper  flyer  on  a  PBX-9404  sample  at  550  n 
the  quadratic  Hugoniot  and  surface  burning  with  v.=0.0794.  The  amount  of  reaction  in  the 
shock  wave  has  been  reduced  but  the  reaction  delay  following  shock  passage  appeals  too 
great. 
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